Universal sub-leading terms in ground state fidelity 
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The study of the (logarithm of the) fidelity i.e., of the overlap amplitude, between ground states of Hamiltoni- 
ans corresponding to different coupling constants, provides a valuable insight on critical phenomena. When the 
parameters are infinitesimally close, it is known that the leading term behaves as 0{L°') (L system size) where 
a is equal to the spatial dimension d for gapped systems, and otherwise depends on the critical exponents. Here 
we show that when parameters are changed along a critical manifold, a sub-leading 0(1) term can appear. This 
term, somewhat similar to the topological entanglement entropy, depends only on the system's universality class 
and encodes non-trivial information about the topology of the system. We relate it to universal g factors and 
partition functions of (boundary) conformal field theory in d = 1 and d = 2 dimensions. Numerical checks are 
presented on the simple example of the XXZ chain. 
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Introduction. — Let |^(A)) denote the ground state (GS) of 
a system with hamiltonian H{X) depending on a set of pa- 
rameters A. We define the ground-state fidelity associated to 
the pair of parameter points A and A' as follows 



,F(A,A'):=|(*(A)|vI/(A'))| 



(1) 



This quantity might provide valuable novel insight for sys- 
tems exhibiting quantum phase transitions IT] |2] |3], in par- 
ticular when there are no obvious local order parameters, but 
some sort of topological order [4|. The strategy advocated 
in Refs Q and |l6| is differential geometric in nature. The 
parameters A and A' are chosen infinitesimally close to each 
other and one focusses on the leading term, the fidelity met- 
ric or susceptibility xl, in the expansion of (fTh as a function 
of (5A := A - A': J^ ~ 1 - SX^XL^f^- Critical lines can 
be identified as singular points of the fidelity metric (in the 
thermodynamical limit) |J5) or by its finite-size scaling |6|. In 
particular in [6 1 it has been shown that the leading finite-size 
term in the fidelity metric is always extensive for gapped sys- 
tems whereas at the critical points its singular part obeys the 
scaling xl/L'^ ^ jjiz+d-2Ax ^ where z is the dynamical ex- 
ponent and A\ the scaling dimension of the operator coupled 
with A. For sufficiently relevant interactions one sees that the 
fidelity metric can display a super-extensive behavior that in 
turns is responsible for the fidelity drops observed at the quan- 
tum phase transition (QPT). On the other hand for marginal 
perturbations Ax ^ d + z i.e., when one is moving along a 
manifold of critical points, the above scaling formula does not 
provide a definite prediction as both log-like terms and 0(1) 
might appear. Accordingly moving along a line of gapless 
point may not give rise to a detectable fidelity drop 16] |3 • 

In this paper we shall demonstrate that the finite-size ex- 
pansion of the GS fidelity (fTl), when A and A' are critical, 
may feature sub-leading terms of order one that depend only 
on the universality class of the considered model and encode 
non-trivial information about the system topology. This goal 



will be achieved by establishing connections to (boundary) 
conformal field theories (BCFT) fSl, I^J. We shall discuss the 
1h-1 free-boson case with the support of exact diagonaliza- 
tion results for the critical XXZ chain, and the 2h-1 quantum 
eight- vertex model |10|. Finally extensions to the case where 
one of the parameters corresponds to a gapped phase will be 
discussed, and potential connections with entanglement mea- 
sures proposed. 

Fidelity and critical theories with boundary. — We would 
like now to establish, on general grounds, a connection be- 
tween the GS fidelity (fTJ and the partition function of a clas- 
sical statistical mechanics system with a boundary interface 
between regions with different coupling A and A'. This can be 
qualitatively understood in terms of the usual correspondence 
between wave functions and path integrals. For the sake of 
concreteness we will now consider space to be compactified 
on a d dimensional hypercube of linear size L (with say pe- 
riodic boundary conditions) and, since we are interested in 
GSs, imaginary time to be infinitely extended along the x^- 
axis. The scalar product of GSs then becomes the (properly) 
normalized partition function for a theory on an infinitely long 
cylinder split into two regions with different couplings A and 
A'. More quantitatively, one can prove the following 



J^(A, A') = lim 



^(A,AO 

Vzwzix') 



(2) 



Here Z{X) is the partition function for the homogeneous sys- 
tem of size 2Lt- X L and Z{X, A') is the partition function in 
the same system with one interface. 

To convince oneself of the validity of Eq. Q one can pro- 
ceed in different ways. Let us assume, for example that the 
underlying classical 2D statistical model, with size L^ x L 
can be described by a transfer matrix T{X). For L^ -^ cx3, 
one has |*(A)) = r(A)^- |$) /^Z(A) with |$) not orthog- 
onal to the ground state and Z{X) the partition function of a 
homogeneous system of size 2L,- x L and boundary condi- 



tions which depend on the quantum model and on |<I>). Then 
(J2]i follows immediately, by introducing the spectral resolu- 
tions of the T's and taking the imaginary time limit Lr ^ oo 
to enforce the projection onto the ground state. Then Z{\, A') 
is the partition function for the system with one interface and 
Z(A) = Z(A,A). 

The sort of inhomogeneous system we have in mind is often 
better seen as a system with a boundary. This is easily done by 
folding: instead of having fields on both sides of the interface 
(where the scalar product is evaluated) one can consider fields 
only on the left side with coordinate Xr < and fold the 
fields living on Xr > into the left domain by introducing 
new species. The problem then becomes a boundary one for a 
theory with double the number of species, and some boundary 
condition (BC) at Xr — 0. 

Boundaries and impurities. — Let us for the moment fo- 
cus on one-dimensional quantum systems d — 1. By again 
using the standard mapping to 2D classical system we have 
In Z{X, A') — InzL — LLrf where / is a non-universal bulk 
term and In zl is a term associated with the boundary itself. 
One can now go back to a d = 1 quantum point of view but 
this time with space along the Xr axis, and L interpreted as 
the inverse temperature f3. One can write the free energy as- 
sociated with the boundary as Lft, := —Ihzl = Lu — s. In 
the critical case for L ^ oo the latter term term gives rise to a 
degeneracy 0(1) factor g = e^, which is, by scaling, indepen- 
dent of L Jl 1 1. This boundary degeneracy - or equivalently s, 
often referred to as the boundary entropy - has played a ma- 
jor role in the analysis of boundary conformal field theories 
(BCFTs). It has been proven in particular that it is universal, 
and thus depends only on the universality class of the critical 
theory, and the type of conformal boundary condition lfT2l : 
for instance, for the Ising universality class with free bound- 
ary conditions, g — 1 , while for fixed boundary conditions 
g — -j^. Note that the issue of scalar product of ground states 
occurred in this context very early on, through considerations 
of the Anderson orthogonality catastrophe lil3J . 

Fidelity andBCFT. — We consider first the archetypal prob- 
lem of a two dimensional free boson with two different val- 
ues of the coupling constants. We write the action as 5* = 
ELi T /d.(^m'^i)' where i?i,2 = R± x [0,L]. The only 
condition we put at the 'interface' x^- = is that the fields 
are continuous (this corresponds to taking the scalar product 
of wave functions). 

We first recall that for a single free boson with coupling 
A, compactified on a circle ip = ip + 2tt, the g factors are 
gjy = 2-i/2(7rA)-i/'' and gN = (TrA)!/^, for Dirichlet (D) 
and Neumann (N) boundary conditions respectively. Assum- 
ing now both bosons compactified on a circle of circumfer- 
ence 2tt, and using the equations of motion to fold the system 
in half gives rise to an equivalent problem of two orthogonal 
species of bosons with the same compactification radius, one 
seeing N boundary conditions with coupling Xn = Ai + A2, 
the other seeing D boundary conditions with coupling A^i — 



7T^^^^. The total g factor is thus 



.9 = .gD(An)gAr(AAr) = 



1 /A 



\/2 VA 



\N 



1/4 



/ Ai + A2 

2v A1A2 



(3) 



Of course we recover that <? = 1 when Ai = A2. Moreover, 
because one field sees D and the other N, it is clear that in fact 
the final result does not depend on the compactification radius 
(and is homogeneous in As). 

It is instructive to recover this result via a direct compu- 
tation (see also fTTi). First note that, in this non-interacting 
case, one can show that Z(Ai,A2) = Z(^J^^). Second, 
since all modes contribute identically to the ratio of partition 
functions (one is simply dealing with Gaussians), the full fi- 

dehty is i^(Ai, A2) = Yik^o V ^7+^ where the product is 
taken over the Brillouin zone. This means L ~ I modes, -the 



zero mode is missing- and thus we have F 



-fL 



g, with g 



the same as (J3]l of course. While in this calculation /, the bulk 
term, is identical to In g, we emphasize that, unlike g, f is not 
universal and depends in general on the details of the model. 
It is interesting to check our prediction against some quick 
numerical calculations. We thus consider XX Z spin chains 
defined on a circle of length L with anisotropy A. Going 
over the standard fermionization and bosonization steps |T5l 
and matching the results with Bethe-Ansatz, one finds that the 
continuum limit corresponds to 



1 r Ai 1 

A = — - \tt — arccosA = . 



(4) 



Here we used the conventions where the spin uf of the spin 
chain is described by dxif, and K is an alternative cou- 
pling constant often used in the condensed matter litera- 
ture. In Fig. [T] we report the results obtained for the lattice 
XX Z model together with the theoretical predictions based 
on BCFT Eq. ([3]) and Eq. (J4]l; the agreement is very good. We 
note that the g factor does depend on boundary conditions. 
For instance it is possible, by breaking the 0(2) symmetry of 
the XXZ chain, to induce antiperiodic conditions on the fields 
(/) in the x direction; a quick calculation shows then that the 
term 0{\) in the fidelity disappears, i.e. 5 = 1 in this case, 
again in agreement with our numerics. 

This kind of calculation admits many variants. Instead of 
having both hamiltonians involved in the fidelity critical, we 
can decide to have only one. In this case, the massive side 
induces a conformal boundary condition on the critical side 
in the calculation of Z(A, A'), and the term of 0(1) in the fi- 
delity is given by the corresponding g factor We can simulate 
this situation by turning again to the XXZ model but this time 
choosing one of the As to be much greater than one. In this 
case, the massive side is in the ordered phase, corresponding 
to two possible ground states, described in terms of spins as 
af — (—1)* and af — (— l)*"*"^ respectively. Each of these 
ground states induces Dirichlet boundary conditions for the 
field ip on the massless side. For each of these Dirichlet cases 
we have gjj — K^^^. Meanwhile, the massive side is a super- 
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FIG. 1 : (color online). Universal g-factor in the fidelity for the XX Z 
model together with BCFT predictions. The fidelity is computed be- 
tween ground states with different anisotropies Ai,2 in the critical re- 
gion (I Ai,2| < 1), and we fixed Ai = 0.20. Plus sign are extrapola- 
tions of data obtained with Lanczos diagonalization on very small lat- 
tices (length L < 22) and periodic boundary conditions. These data 
agree perfectly with BCFT predictions Eq. (|3l together with Eq. J4l. 
Instead, toroidal BC given by cr|^, j^ = e^ af, af^^i — ~af in- 
duce antiperiodic BC on the field and no g factor (i.e. g — 1). Note 
that such BC (they belong to conjugacy class (IV) of Ref. 1 16|) break 
conservation of total magnetization. In the inset the fidelity is com- 
puted when one ground state is critical and the other is deep the mas- 
sive (Neel) phase Ai ^ 1. Solid line give the BCFT prediction. The 
small discrepancy around A2 ~ 1 is due to finite size effects which 
are more pronounced near the Kosterlitz-Thouless point A = 1. 



position of the two orthogonal ground states with equal coeffi- 
cients 4= so we get in the end g = 2x-y=x K'^/^ = y/2K^/^. 
Again these predictions are well confirmed by finite size Lanc- 
zos calculation (see inset of Fig.fT]). 

Terms of order one in the 2+1 case: the quantum eight- 
vertex model — We turn to consider 0(1) terms in the GS fi- 
delity of 2 + 1 models whose quantum critical points have 
dynamical critical exponent z = 2. For these models at criti- 
cality, ground state functionals are conformal invariant in the 
2d physical space, and equal time coiTelators coincide with 
correlations in a 2d CFT. We will now show that the fidelity 
involves universal terms of 0(1) in this case as well, and that 
this time they are related to partition functions of CFTs on 
Riemann surfaces. 

To make things concrete, let us specialize to the 2+1 ana- 
log of the free boson - the quantum Lifschitz model - for 
which a convenient lattice realization is provided by the quan- 
tum vertex model ifTOl . The Hilbert space of this model is 
spanned by an orthonormal basis {|C)} in one-to-one corre- 
spondence with the configurations of the classical eight- vertex 
model. The hamiltonian has the form H = ^^ Qi, with Qi 
positive operators, chosen such that H annihilates the follow- 
ing state: \^{c^)) = J2{c} c^^^'^^C)/^Z2d{c^), where we 
have chosen, for simplicity, a = b=l,d = Oso the only 
remaining parameter is c, which is the equivalent here of A 
in the previous sections. The normalization factor is given 
by the partition function of the classical eight-vertex model: 
Z2d{c^) = S|c> c^"'=^''\ where hc{C) is the number op- 
erators for the ctype vertices, for the configuration C. The 




FIG. 2: (color online). Universal g-factor in the fidelity of the quan- 
tum eight-vertex model with periodic BC when all the theories are 
in the disordered region, i.e. c* < 4, (c')* < 4, {cc')^ < 4. We 
fixed Li — L2. The g factor is smooth at the border of the region 
c, c' ^ and c, c' -^ \/2- 



ground-State phase diagram for the quantum model is iden- 
tical to the classical one, but given in terms of c^. The scalar 
product of ground states is given by 



(*!*') 



Z2d{cc') 



V^2d(c2)Z2d((c') 



(5) 



The models live in the 2d equal time slice, and as usual we 
will take the dimensions Li,L2 of this slice to infinity. Now 
consider the case where the weights obey c^ < 2, (c')^ < 
2, \cc'\ < 2. In that case we are dealing with paitition func- 
tions of two dimensional critical vertex models, which are de- 
scribed in the continuum limit by Euclidian free bosons in 2d 
L8J . With periodic boundary conditions for instance, these 
partition functions behave as Z2D = g~^^^^'^Zcft{Li/ L2) 
where Zqft is the modular invariant partition function of the 
conformal invariant field theory. 

The important point is that the scalar product (Bl will 
have a term behaving like an exponential of the area 
exp {-L,L2 [f{cc') - /(c2)/2 - /((c')2)/2] } and 
of order one, Z{A)Z-'^^^{X)Z-'^^^{\') where A 
mally the same coupling constant as before (Hli 
is the coupling associated with the product cc': 
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is for- 
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- 1 = — cos27r^A, 
= — cos27r^A. The conformal partition function 

reads Z{X) = {ri{q)ri{q))^^T{X) where: q 

- exp(— 27rLi/L2) parametrizes the torus, 1(A) 



Eoo 
n,m—- 



A 



H-. 



= +mV27rA)^ -|( ^ji^ -mV27rA)^ 



and 77(g) — 



9^ ni^i(l ~ 9") ID- While the prefactor and the 77 terms 
disappear in the ratio, the instanton sums T remain, leading 
to an rather complicated expression X(cc')/-\/X(c^)X((c')2) 
for the term O ( 1 ) . An example of the behavior of this term is 
given in figure l2] 

Interestingly, being a combination of conformal partition 
functions, the term of 0(1) depends heavily on the topol- 
ogy and boundary conditions of the base space. One can for 
instance imagine defining the 2d quantum models on higher 
genus Riemann surfaces II17I . or on surfaces with boundaries 
and curvature. To give a very simple example, the logarithm 



of the free-boson partition function on a rectangle of size 
Li , L2 with free boundary conditions (either D on all sides 
or N on all sides) is given by \nZ2D = aLiL2 + h{Li + 
L2) + J In L2 — I In ['7(9)] where a, h are non universal terms. 
The logarithmic term meanwhile is universal, and comes from 
the general formula for the free energy of a critical region A 
of linear size L2 oc Li, Euler characteristic x ^nd a boundary 
with a discrete set of singularities. Note that all dependence 
on the coupling constant A has disappeared! It follows that if 
we were to calculate the scalar product of ground states in this 
situation, there would be just no term of order one. 

Connections with quantum entanglement. — Before con- 
cluding we would like briefly to comment about the possible 
connections between the fidelity approach pursued in this pa- 
per and quantum entanglement. First let us notice that BCFT 
arguments have been used in the calculations by Kitaev and 
Preskill to motivate their expression for the topological entan- 
glement entropy (TEE) 1 18 1. Their derivation shows the TEE 
is a 0(1) sub-leading universal term that is strictly analogous 
to those for the fidelity in this paper This is even more man- 
ifest if one expresses the degeneracy g factors in terms of the 
modular S matrix of the CFT 1 1 1 1 and compares it with the 
TEE Stopo = log l^*"!. Moreover there is a striking similarity 
between our formulas for the log of the fidelity in the previous 
section and formulas in 1 19| for the entanglement entropy at 
conformal quantum critical points. Both involve logarithms 
of conformal partition functions, and it is clear that, by taking 
the ground states of the quantum vertex model with different 
couplings in different regions, one could obtain entanglement 
entropy through a term of 0(1) in the fidelity. How general 
and useful this observation might be is an open question. 

Conclusions. — Using BCFT techniques, we have shown 
that the fidelity between critical states contains a term of or- 
der 0(1) which depends only on the universality class and on 
the topology of the base space. As such it bears similarity to 
the topological entanglement entropy or the central charge ap- 
pearing in the expansion of the ground state energy. The use 
of methods of CFT in information theory should go much be- 
yond the consideration of these terms of 0(1). For example, 
the same techniques can be used to extract information about 
the Loschmidt echo Ii20il . 
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